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In a recent paper [1] we have shown that the General Relativistic (GR) 
dynamics of a self-gravitating perfect fluid is greatly simplified under three 
assumptions: i) the fluid is coUisionless (i.e. with zero pressure, p), ii) it has zero 
initial vorticity, ojab [2] and iii) the so-called "magnetic" part of the Weyl tensor, 
Haft, is zero. The former two conditions are wide enough to allow for many 
cosmological cases, such as the evolution of dark matter adiabatic perturbations 
generated during inflation. The third assumption is more problematic. In linear 
theory Hab only contains vector and tensor modes (e.g. Ref.[3]): if the vorticity 
vanishes no vector modes are present and B.ab only contains gravitational waves. 
Beyond linear theory the meaning of is less straightforward. It is reasonable 
to assume that Hab = forbids at least the occurrence of gravitational waves. 
This is particularly clear in the present context, where, thanks to the absence of 
pressure gradients, the motion is geodesic and, if Hab also vanishes, no spatial 
gradients appear in the evolution equations (apart from those contained in 
convective time derivatives, which can be dropped by going to a comoving 
frame): it is hard to think of any actual wave propagation with no spatial 
derivatives appearing in the fluid and gravitational evolution equations. 

Following Ellis [4] we describe the dynamics directly in terms of observ- 
able fluid and geometric quantities: the mass density g, the expansion scalar 
and three traceless, flow-orthogonal and symmetric tensors, the shear, cr'*^, the 
so-called "electric" part of the Weyl tensor, E'*^, describing tidal interactions of 
the fluid element with the surrounding matter, and its magnetic part H'*^. As 
noted in Ref.[l], if the magnetic component is switched off, all the equations for 
the GR dynamics take a strictly local form: each element evolves independently 
of the others. Only at the initial time Cauchy data must be consistently given 
on a spatial hypersurface. The subsequent evolution can be entirely followed 
in Lagrangian form until caustic formation, when the one-to-one mapping be- 
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tween fluid elements and space points is lost. We call such a system a silent 
universe, in that no information can be exchanged among different fluid ele- 
ments: this is due to the causal nature of GR, where signal exchange can only 
occur dynamically via gravitational radiation and, in the case of fluids with 
non-zero pressure, also via sound waves, but none of these wave modes is al- 
lowed when p = = 0. Because of the advantages of a purely local treatment, 
this method [1] has recently attracted some attention. In particular, Croudace 
et al. [5] have shown the connection of the GR pancake solution [1] with the 
Szekeres metric [6]; Bertschinger and Jain [7] have performed a detailed study 
of the Lagrangian dynamics of fluid elements. 

However, the condition Hab = cannot be taken as an exact constraint 
for the general cosmological case. It has been shown [8] that the only solutions 
of Einstein equations, with p — Uab — i^ab — are either of Petrov type I, or 
conformally flat, or homogeneous and anisotropic of Bianchi type I, or locally 
axisymmetric (i.e. with two degenerate shear eigenvalues) and described by a 
Szekeres line-element [6]. All of these cases require some restrictions on the 
initial data: the exact conditions above are not suitable to study cosmological 
structure formation. However, requiring p = u>ab = and H^b ~ appears more 
feasible. A small Kab is in fact compatible with arbitrary departures from local 
axisymmetry of fluid elements. This is shown by the behaviour of perturba- 
tions around Robertson- Walker (RW): whatever initially scalar perturbations 
are given, Kab vanishes at first order, but not beyond. A small value of B.ab 
allows arbitrary ratios among the shear eigenvalues, provided the initial pertur- 
bations are small. For general initial shapes of the fluid elements the system will 
radiate gravitationally during non-linear evolution. However, fully GR numer- 
ical computations [9] have shown that only a negligible fraction (less than 1%) 
of the total energy is carried away in the form of gravitational radiation, during 
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the non-linear collapse of coUisionless ellipsoids. In spite of these facts, as our 
calculations below demonstrate, a non-zero H^b allows for the influence of the 
surrounding matter on the evolution of fluid elements. Although this signal 
travels at finite speed, for perturbations on scales much smaller than the hori- 
zon it effectively appears as an instantaneous Newtonian feature. One might 
wonder whether during the late phases of collapse, when local axisymmetry is 
expected to be established, the environmental influence on the evolving fluid 
element can be neglected and the Hab = condition restored. 

General relativistic dynamics - To describe our system we start from the 
equations of Ref.[4]. We always work in the comoving synchronous gauge 
ds^ = —dt^ + a^{t)ha/3dq°'dq^, where a = At"^/^, as for a flat, matter- 
dominated RW model (our "background" solution). For computational con- 
venience we introduce suitably rescaled quantities: a scaled density fluctuation 
A = [QTrGt^Q — l)/a, a peculiar expansion scalar d = (3t/2a)(0 — 2/t), a trace- 
less shear tensor s"^ = (3t/2a)(j"^ and a traceless tidal tensor e*^^ = (3t^/2a)E"^. 
These quantities can be grouped in two space-like tensors: the velocity gradient 
tensor •d°'p = s'*^ + ^S°'f3'&, related to the covariant derivatives of the peculiar 
velocity field; the peculiar gravitational field tensor A"^ = e"^-|- lA^'^. We also 
scale the magnetic tensor as = (3t^/2a)H"^. The dynamical equations for 
the fluid and the gravitational fleld are 

A a 

n p 



= --(^«^ + A-^)-^"^^-^, (1) 

+ 5«^(A^ - A^,^^^) + —hpr^iv^-^'n'^^-s +r'''n\,s ) , (2) 

- T^hf,r,{fj^'^'A-^;s +r'''^\,5 ) , (3) 



where the dot denotes partial differentiation with respect to the scale factor a 
and fj'^'^^ is the Levi-Civita tensor relative to the metric ha/3: fj"'^'^ = h~^/'^e°''^^ , 
with £^^^ = 1. The metric tensor evolves according to ^h'^'^h^p = ■i?"^. 
The above tensors have to satisfy the constraints [4] 

^/;^=^,a, (4) 

A£;p = A,« -^ha^hp,fi'^^^^\'Hf'^ , (5) 
= —ha^hf^.r'"' ^\A^-y , (6) 

= ^hAv''^'K-^5 +r^'^^^;s ) • (7) 

All these are fulfilled at the linear level [3] by growing-mode scalar initial condi- 
tions [1]: A"^(ao) = — 'i?"^(ao) = <fQ"i3, where the scalar (fo, an arbitrary func- 
tion of the space coordinates q°', is the initial peculiar gravitational potential, re- 
lated to Bardeen's gauge-invariant $h [10] by (fo = — (3/2A^)$//. These initial 
conditions correspond to the "seed" metric hap = 5a/3(l — ^^^V'o) — 2a<^0)a/3) 
and imply vanishing initial 'H°'p (the constant mode, oc A^(fQ <^ 1, can be 
neglected in practice, compared to the growing mode cx ai/Jca/?)- 

The Lagrangian dynamics is determined by Eqs.(l), (2) and (3) plus the 
initial data. One obtains a local Eulerian description of the fluid [1], using 
the "generalized Hubble law" [4]. We have = where a$,°' is the in- 

finitesimal spatial displacement of neighbouring elements. The matrix connect- 
ing the Eulerian coordinates x'^ with the Lagrangian ones is the Jacobian 
J"^ = dx'^/dqf^ = + P"^, where X>"^ is the (symmetric) deformation tensor. 
Taking = dx"^ = J^jj^^Q^ where ^^^^ = dq^ represent the initial (i.e. La- 
grangian) infinitesimal displacements, one gets 1)°'^ = i?"^ -|- ■&°'^T>'^^, formally 
solved by V°'p{a) = exp dai!}'^p{a) — 5"p. Once the Jacobian is known one 
gets the metric as hajs = h^5{aQ)J\J^p. As shown in Refs.[8,l], if = 0, the 
tensors i?"^, A"^, hap commute and they can be diagonalized simultaneously. 
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In such a case, Eqs.(l) and (2) can be reduced to six first order equations for 
the six eigenvalues of t?"^ and A"^^. Along the local principal axes we can set 
haf3 = ^aphp and = S'^p-djs and get ha{a) = /ia(ao)exp2 j^^da{}a(a). In 
the locally axisymmetric case, i.e. when two eigenvalues of ipoi'^p coincide, a 
relation exists with particular Szekeres solutions [6] . 

Newtonian dynamics - The equations which govern the non-linear dynamics of 
a coUisionless fluid in Newtonian Theory (NT) for an expanding universe [11] 
can be written in suitably rescaled form as (e.g. Ref.[12]) 

+ A«,^=-|-(ti« + ^,«) , (8) 
la 

A + .x/3A,^=-i(«^^+A)-A..^^ , (9) 
a 

(10) 

where is the peculiar gravitational potential. Differentiating the Euler equa- 
tion (8), defining the symmetric tensors = u°'.,f3, with u°' = dx'^/da, and 
A"^ = v'j'^/jj ciiid adopting a Lagrangian description, one recovers Eq.(l), while 
the continuity equation (9) coincides with the trace of Eq.(2). It is clear that 
the NT is degenerate, as it provides only one equation to determine the tensor 
A°^: any traceless tensor added to the r.h.s. of Eq.(2), leaves the NT equations 
unchanged. In order to completely determine the evolution of the gravitational 
field tensor A'^p one has to resort to its definition in terms of the potential (p, 
i.e. to a non-local theory. Because of the intrinsic non-locality of NT (the 
Poisson equation (10) is an elliptic, constraint equation) one needs boundary 
conditions to determine the dynamics: contrary to the GR equations, initial 
data are not enough. It is well known (e.g. Ref.[4]) that the lack of evolution 
equations for the traceless part of the gravitational tensor, e"^, implies that the 
NT adds spurious solutions which would be discarded by the full GR system. 

6 



Beyond the Zel'dovich approximation - In order to see the behaviour of the GR 
solutions and evaluate the role of the magnetic term we construct a second order 
Lagrangian perturbation expansion in the amplitude of the fluctuations around 
RW. It will prove useful to define the two quantities jii = (po,'^^ = Ai + A2 + A3 
and fj,2 = |(vo,\ V^o,''^ -'Po,'^s '^o,\ ) = A1A2 + A1A3 + A2A3, where Aq, are the 
local eigenvalues of the symmetric tensor <fQ,'^i^- One immediately obtains the 
traces, = —jii + a{—ji\ + and A = /xi + a{ji\ — ^IJ,2), which coincide 
with those obtained in Lagrangian second order NT [13]. After very lengthy 
calculations we obtain 



= -^o,"^+^(-12//l^O,"^+6//2<5'^;3 + 5^o,"^^0,^^) +X"/3 , (11) 



(here indices are raised by the Kronecker symbol), having kept only grow- 
ing modes. The expressions for A"^ and 7Y"^ will not be reported here 
for shortness. The traceless tensor x"^, representing the contribution due 
to the magnetic part, has zero divergence: X/s'ja = 0. It can be writ- 
ten as a convolution x°'f3i<li — J '^'^'?'5'"^(q')/(|q — q'|,a), of the source 
S'^p = A*2,"/3 -V^(2//i(y?o,"/3 -2<^o,"y <y^'o,^/3 -^"^A*2) with the function /, whose 
Fourier transform f{k) satisfies the equation, 

where a prime denotes differentiation with respect to the conformal time r = 
{3/A)t^/^. The initial conditions are /(tq) = /'(tq) — /"(tq) = 0. Asymptotic 
solutions of Eq.(12), confirmed by a numerical check, are / ~ 2A^t^ /21k^ for 
kr ^ 1 and / fa for kr <^ 1. Performing a second-order expansion 

for the deformation tensor and defining = — a(/Po,"/3 +'^^'0'^; we find 
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with trace ip'^^ = — |/i2- The symmetric tensor 'ip°'p provides the second or- 
der correction to the deformation tensor, whose first order is the kinematical 
Zel'dovich approximation. The metric tensor reads 

haf3 = 5al3 " 2aV3o,a/3 +y ^19<y^'0,a7 ¥^0,^/3 -12//l¥'0,a^ +6//2(^a/3^ + j daXa/3 ■ 

(14) 

We then have dx" = dg" — a(po,'^pdq^ + a^ip'^j^dq^. In NT one would 
write the same formal expression, but the irrotationality condition would lead 
to = '0,'*^, with the potential ip satisfying the second order Poisson equation 
[13] V^V' = ~f A*2) which is consistent with the trace of the GR equation. In 
other words, the NT eigenvalues Ua of 'ip°'fj only need to satisfy the condition 
= — |/i2- In order to get the complete information on the single v^s 
one needs the NT definition of '0"^ as '0;"/3; a non-local information. The 
GR also solve the NT equations, but the reverse is not necessarily true: it 
depends upon the boundary conditions used in solving Poisson's equation. 

Inside the horizon - Suppose that the source, hence (fo,°'/j, has some typi- 
cal scale of variation £, i.e. £ ~ (po^^^f^ /(po,°'p^- If ^ <S r we find ^ 
—(fo,'^l3—0''fio,'^'yVo,^f3+'20'4'"i3- The second order deformation tensor reduces 
to ip'^p = ipj'^fj, while the metric reads ha/s = Sa/3 — 2a(fo,a/3+a^'4'jai3- All these 
expressions coincide with those of second order NT and can be obtained from 
the c — > oo limit of Eq.(12). The scalar ijj carries information on the infiuence 
of the surrounding matter on the dynamics of fluid elements. Note that 5*^/3 
produces a tilt of the principal axes of the first-order deformation tensor, (foj'^p- 

Outside the horizon - When £ :$> r, x'^js ~ (3i^/14a)»S'^, and the contribution 
to d'^p due to the magnetic term becomes negligible. The relevant expressions 
can be obtained from Eqs.(ll), (13) and (14) with 0. Perturbations with 

size greater than the Hubble radius evolve as a separate silent universe: spatial 
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gradients play no role in this case. However, these local GR effects have little 
cosmological implications, since perturbations on super-horizon scales usually 
have very small amplitude, and a linear approximation is sufficient. Neverthe- 
less, there are a number of formal consequences, which is worth mentioning. 
One of these is the absence of 2D solutions. If one eigenvalue of <^0)"/3) e.g. 
A3, vanishes everywhere, the NT, with suitable boundary conditions, implies 
'dz{a) = or ^3(0) = q'3, i.e. no motion along the third axis. This is referred as 
"two-dimensional" (2D) gravitational clustering. As far as the second order de- 
formation tensor is concerned, one would have 1^1 + ^2 = — f A*25 with jj,2 = A1A2, 
and = 0. In the GR case, instead, we find ui = 1^2 = —I's = — f A*2) and 
^3(0^) 7^ for o 7^ ttQ. The motion dynamically impressed along the third axis 
soon becomes of the same order of magnitude as that in the other directions. 
This effect is due to the tide-shear coupling term 5°g('!?A — A'^g-d^^) in the evo- 
lution equation for A"^, which reduces to — 2(U25"^ to lowest order. The only 
case when this coupling disappears is when two A^'s simultaneously vanish, i.e. 
for planar symmetry. Therefore 'i?3(a) = is not an exact solution of the GR 
equations, unless another -da also vanishes. As an example, no axisymmetric 
configurations without motion along the symmetry axis are allowed. 

This discussion leads to the main issue: the general non-linear dynamics 
of fiuid elements. So far, two analytical solutions of our system are known: for 
planar configurations, Ai = A2 = 0, one recovers the Zel'dovich pancake solu- 
tion, as shown in Ref.[l]; for exactly spherical configurations, Ai = A2 = A3, 
the local solution is the well-known top-hat model (e.g. Ref.[ll]). Croudace et 
al. [5] looked for solutions representing attractors among the trajectories of our 
system with zero magnetic tensor. They found that both spherical collapse and 
a perfect pancake are repellers for general initial conditions, and argued that 
the pancake instability is probably due to having disregarded the contribution 
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of 'H°'p. On the other hand, Bertschinger and Jain [7] have shown that the 
instabihty of the pancake solution is caused by the tide-shear couphng in the 
evolution of the tide, which tends to destabilize the pancake solution (for gen- 
eral initial conditions) but would stabilize prolate configurations. For vanishing 
Ti-'^pi a strongly prolate spindle with expansion along its axis is the generical 
outcome of collapse, except for specific initial conditions corresponding to ex- 
actly spherical or planar configurations. Our analysis shows that the dynamical 
effect causing preferential collapse to expanding spindles in the Ti'*^ = case 
is the GR tide-shear coupling in the tide evolution equation. This term is not 
present in NT, although it is compatible with its equations. This is further 
illustrated by the collapse of an infinite homogeneous ellipsoid — > cxo), which 
is described by our equations with zero H"^. As well known, the NT dynamics 
favours the formation of oblate spheroids (e.g. Ref.[14]), pancake-like objects 
with one collapsing axis and the other two tending to a finite size (apart from 
initial conditions corresponding to an initial prolate spheroid). The GR collapse 
(e.g. Ref.[15]) favours the formation of prolate spheroids, collapsing filaments 
with expansion along their symmetry axis. However, as our second order calcu- 
lations show, the evolution of fluid elements as isolated ellipsoids does not apply 
to perturbations on scales smaller than the Hubble radius: here non-local effect 
play a fundamental role. The actual non-linear dynamics would generally re- 
sult from the competition of the local GR tide-shear coupling, causing pancake 
instability, and the non-local environmental influence, carried by the magnetic 
part of the Weyl tensor. During the early deviations from linear evolution, 
as described by Lagrangian second order perturbation theory, the latter effect 
dominates; however, extending this conclusion to the late strongly non-linear 
phases would require further study. 

Finally, as a result of this analysis, we are able to calculate how many 
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gravitational waves are produced within a second order approximation (re- 
member that at this order the magnetic tensor is traceless and transverse, 
so it is related to gravitational radiation [3]). Outside the horizon B.a/3 ~ 
{4:a^/7t){eajd'Po,^p'Po,\+£i3j6'Po,^a'Po,\),'', while inside the horizon Hq^ de- 
cays as 1/a. Only a tiny amount of gravitational waves is produced on sub- 
horizon scales at this level, however their dynamical role is far from being neg- 
ligible! 
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